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ABSTRACT 



1-^ ■ Images of the site of the Type Ic Supernova 2002ap taken before explosion were 

analysed previously by Smartt et al. (2002). We have uncovered new unpublished, 
Q , archival pre-explosion images from the Canada-France-Hawaii Telescope (CFHT) that 

■ are vastly superior in depth and image quality. In this paper we present a further search 

c/3 \ for the progenitor star of this unusual Type Ic supernova. Aligning high-resolution 

C$ . Hubble Space Telescope (HST) observations of the supernova itself with the archival 

CFHT images allowed us to pinpoint the location of the progenitor site on the ground 
\ based observations. We find that a source visible in the B and R band pre-explosion 

J> ■ images close to the position of the SN is (1) not coincident with the SN position 

' within the uncertainties of our relative astrometry, and (2) is still visible ~4.7 yrs 

<~) post explosion in late-time observations taken with the William Herschel Telescope. 

IO ' We therefore conclude that it is not the progenitor of SN 2002ap. We derived abso- 

lute limiting magnitudes for the progenitor of Mb ^ -4.2 ± 0.5 and Mr ^ -5.1 ± 0.5. 
\Q . These are the deepest limits yet placed on a Type Ic supernova progenitor. We rule 

f^i ' out all massive stars with initial masses greater than 7-8M Q (the lower mass limit 

\ for stars to undergo core collapse) that have not evolved to become Wolf-Rayet stars. 

This is consistent with the prediction that Type Ic supernovae should result from the 
explosions of Wolf-Rayet stars. Comparing our luminosity limits with stellar models of 
single stars at appropriate metallicity (Z=0.008) and with standard mass loss rates, we 
find no model that produces a Wolf-Rayet star of low enough mass and luminosity to 
be classed as a viable progenitor. Models with twice the standard mass loss rates pro- 
vide possible single star progenitors but all are initially more massive than 3O-4OM . 
We conclude that any single star progenitor must have experienced at least twice the 
standard mass loss rates, been initially more massive than 30 — 40-Mq and exploded 
as a W-R star of final mass 10-12M Q . Alternatively a progenitor star of lower initial 
mass may have evolved in an interacting binary system. Mazzali et al. (2002) propose 
such a binary scenario for the progenitor of SN 2002ap in which a star of initial mass 
15-2OM is stripped by its binary companion, becoming a 5M Q Wolf-Rayet star prior 
to explosion. We constrain any possible binary companion to a main sequence star 
of ^ 2OM0, a neutron star or a black hole. By combining the pre-explosion limits 
with the ejecta mass estimates and constraints from X-ray and radio observations we 
conclude that any binary interaction most likely occurred as Case B mass transfer, 
either with or without a subsequent common envelope evolution phase. 

Key words: stars: evolution - supernovae: general - supernovae: individual: 
SN 2002ap - galaxies: individual: M 74 
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1 INTRODUCTION 

Stellar evolution theory predicts that stars initially more 
massive than 7-8Mq should explode as Core Collapse Super- 
novae (CCSNe) at the end of their nuclear burning lives (e.g. 
Heger et al. 2003; Eldridge & Tout 2004). The observation 
of progenitor stars of CCSNe in archival pre-explosion im- 
ages, provides us with a direct means with which to test the 
theoretical predictions. Until quite recently the only directly 
observed progenitor stars of CCSNe were those of the Type 
II-peculiar SN 1987A in the Large Magellanic Cloud (LMC) 
(Walborn et al. 1987) and the Type lib SN 1993J in M 81 
(Aldering et al 1994). The advent of large data archives, per- 
haps most importantly that of the Hubble Space Telescope 
(HST), has allowed several groups around the world to ex- 
tend the progenitor search to larger distances (e.g. Smartt 
et al. 2004; Li et al. 2005; Gal- Yam et al. 2007). To date 
a total of nine progenitor stars of CCSNe have been re- 
ported, all of which are for hydrogen-rich Type II SNe, and 
with the seven most recent additions specifically for Type 
II-P SNe. There has been as yet no direct detection of a 
hydrogen-deficient Type lb or a hydrogen and helium defi- 
cient Type Ic supernova progenitor, although the precursor 
of the peculiar SN 2006jc was observed some two years prior 
to explosion during an LBV-like outburst (Pastorello et al. 
2007). Upper limits have been determined for the luminosity 
of several Type Ib/c progenitors, of which the most restric- 
tive limits set thus far are for the progenitor of the Type Ic 
SN 2004gt (Maund, Smartt & Schweizer 2005; Gal- Yam et 
al. 2005). In terms of absolute magnitude, the pre-explosion 
HST observations probed down to My — -5.3. This allowed 
both groups to determine that the progenitor was almost 
certainly a Wolf-Rayet star, a massive star that has lost its 
hydrogen-rich envelope. Such stars are predicted to be the 
progenitors of Type Ib/c SNe. 

Another Type Ic SN with archival pre-explosion obser- 
vations is SN 2002ap. It was discovered by Y. Hirose on 2002 
January 29.4 UT at magnitude V=14.54 in the spiral galaxy 
M 74 (Nakano et al. 2002), at a distance of approximately 9.3 
Mpc (Hendry et al. 2005). The object was quickly revealed 
as a Type Ic SN with broad spectral features, indicative of 
very high velocities of the SN ejecta (e.g. Meikle et al. 2002; 
Kinugasa et al. 2002; Gal- Yam et al. 2002; Mazzali et al. 
2002; Foley et al. 2003). Although its spectra appeared sim- 
ilar to that of the peculiar SN 1998bw, it was less luminous, 
reaching a peak magnitude of My — -17.5, about 1.7 mag 
fainter than SN 1998bw. Furthermore, unlike SN 1998bw no 
gamma ray burst (GRB) was detected coincident with the 
position of SN 2002ap (Hurley et al. 2002a). UBVRIKHa 
pre-explosion observations of the site of SN 2002ap in M 74, 
taken with the KPNO 0.9 m, the 2.5 m Isaac Newton Tele- 
scope (INT) and the Bok 2.3 m telescope, were used by 
Smartt et al. (2002) to derive upper limits to the absolute 
brightness of the progenitor star in the various filters. The 
B and V band observations were the deepest, reaching ab- 
solute magnitude limits as faint as Mb = -6.3 and Mv = 
-6.6. 

At the time it was believed that these were the 
best quality images available of the pre-explosion site of 
SN 2002ap. However we have discovered additional pre- 
explosion images of the SN position in the archive of the 
Canada France Hawaii Telescope (CFHT), which are vastly 



superior in depth and resolution to those originally pre- 
sented in Smartt et al. (2002) . In fact we shall show that they 
are the deepest pre-explosion images of any nearby Type 
Ib/c SN to date. The new archive at the Canadian Astron- 
omy Data Centre (CADCj3 re turns these images. This paper 
presents the study of these very deep, high quality ground 
based observations of the pre-explosion site of SN 2002ap, 
in which we attempt to detect the progenitor star and place 
rigorous constraints on its properties. 



2 OBSERVATIONS 

2.1 Pre-explosion observations 

The galaxy M74 was observed with the 3.6-m Canada France 
Hawaii Telescope in October 1999, and to our knowledge the 
observations have not been published. The site of SN 2002ap 
was imaged using CFH12K, a 12 CCD mosaic camera, which 
is now decommissioned. The pixel size of this instrument 
was 0'.'2 and further details can be found on the CFHT web- 
sitfi The B and R observations were particularly deep, con- 
sisting of individual, dithered 600-second exposures which 
combined to result in 1 hour 20 mins and 1 hour exposure 
times respectively. The V and Ha observations were much 
shallower with the individual frames combined giving to- 
tal exposure times of just 600 and 900-seconds respectively. 
A summary of these observations can be found in Table [TJ 
The images were downloaded from the CADC archive, and 
most had been processed successfully through the ELIXIR 
system. However, some of the R band images generated by 
ELIXIR were corrupted, so we used the master flat from the 
archive and processed the frames manually using standard 
techniques in IRAF. 

There were 8 x 600s exposures from 1999 October 9th 
in B and the combined image quality was 0'.'75. A further 
8x600s set of exposures were taken on 1999 December 31st, 
but seven of these had an image quality of ~l'/2 (the re- 
maining image had poor seeing of l'.'9). Combining the l'/2 
stacked image with the better seeing (f.'75 images did not 
increase the sensitivity, so we chose to proceed with analysis 
of the 1999 October 9th images only. The 6x600s R band 
images taken on 1999 October 6th were all of similar resolu- 
tion (Cf.'7), and the shorter 5 x 120s V stack had a resultant 
quality of l'.'2. There were 5x300s Ha exposures taken on 
1999 Jul 13, each with image quality around (f.'9. The ob- 
servations were made towards the end of the night which re- 
sulted in the final two exposures having significantly higher 
background counts. These two exposures were excluded from 
the final stacked image since they would only serve to add 
excess noise. 

To determine zeropoints and colour corrections for each 
filter, magnitudes of standard stars identified from the cata- 
logue of Henden et al. (2002) were measured using the aper- 
ture photometry task within the iraf daophot package. 
Due to the depth of the observations many of the standard 
stars in the field were saturated and it was only possible 
to find five useable objects common in BVR. These stars 



1 http://cadcwww.dao.nrc.ca/cfht/ 

2 http:/ /www. cfht.hawaii.edu/Instruments/Imaging/CFH12K/ 
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Table 1. Summary of observational data. 



Date of Observation 


Telescope + 


Filter 


Exposure Time 


Observer 




Instrument 




(sec) 




Pre-explosion observations 










1999 Jul 13 


CFHT+CFH12K 


Ha 


900 


Cuillandre, Hainaut, McDonald 


1999 Oct 06 


CFHT+CFH12K 


R 


3600 


Astier, Fabbro, Pain 


1999 Oct 09 


CFHT+CFH12K 


B 


4800 


McDonald, Cuillandre, Veillet 


1999 Oct 09 


CFHT+CFH12K 


V 


600 


McDonald, Cuillandre, Veillet 


Post-explosion observations 










2003 Jan 10 


HST+ACS/HRC 


F475W 


120 


Kirshncr 


2003 Jan 10 


HST+ACS/HRC 


F625W 


120 


Kirshncr 


2003 Jan 10 


HST+ACS/HRC 


F658N 


800 


Kirshner 


2003 Jan 10 


HST+ACS/HRC 


F814W 


180 


Kirshner 


2004 Jul 06 


HST+ACS/HRC 


F555W 


480 


Filippenko 


2004 Jul 06 


HST+ACS/HRC 


F814W 


720 


Filippenko 


2004 Aug 31 


HST+ACS/HRC 


F435W 


840 


Filippenko 


2004 Aug 31 


HST+ACS/HRC 


F625W 


360 


Filippenko 


2006 Aug 25 


WHT+AUX 


B 


4800 


Benn, Skillen 


2006 Aug 27 


WHT+AUX 


R 


3600 


Benn, Skillen 



CFHT = 3.6-m Canada France Hawaii Telescope, Mauna Kea, Hawaii 
HST+ACS/HRC = Hubble Space Telescope + Advanced Camera for Surveys/High Resolution Camera 
WHT = 4.2-m William Herschel Telescope, La Palma, Canary Islands 



Table 2. Standard star magnitudes taken from Henden et al. 
(2002). 



Star Henden catalogue B V R 

identity 

A 118 20.827 19.286 18.257 

B 105 20.231 19.756 19.523 

C 91 20.804 19.286 18.234 

D 87 21.333 19.827 18.696 

E 85 20.739 20.536 20.326 



Table 3. CFHT zeropoints and colour corrections. 



Filter Colour Correction Zeropoint 
(Colour) 

B 0.04(B-R) 26.10 

R 0.02(B-R) 26.29 



nearby SNe. Details of these observations are given in Table 
[1] All HST observations were downloaded from the Space 
Telescope Science Institute (STScI) archive via the on-the- 
fly re-calibration (OTFR) pipeline. Those from 2004 were 
found to be well aligned despite being taken on two dif- 
ferent epochs almost two months apart. Identification of 
SN 2002ap in these late-time images was aided by compar- 
ison with those from January 2003 in which the SN was 
significantly brighter. In all cases the pointing of the obser- 
vations was such that SN 2002ap was imaged very close to 
the centre of the HRC chip. 

Photometry was carried out using the ACS module 

© 2007 RAS, MNRAS 000, 000-000 




Figure 1. Standard stars in CFHT pre-explosion observations. 
The approximate location of the site of SN 2002ap is also indi- 
cated. All other images in this paper have the same orienation as 
this figure. 



are shown labelled in Figure \T\ and the corresponding cata- 
logue magnitudes are presented in Table [2] The calculated 
zeropoints are shown in Table [3] 



2.2 Post-explosion observations 

Observations of SN 2002ap were taken using the High Reso- 
lution Channel (HRC) (pixel scale = / .'025) of the Advanced 
Camera for Surveys (ACS) on board HST during January 
2003, and July and August 2004. These formed part of the 
programmes G09114 (PI: R. Kirshner) and SNAP10272 (PI: 
A. Filippenko), designed to obtain late-time photometry of 
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of the PSF-fitting photometry package DOLPHOl|f|, a modi- 
fied version of HSTphot (Dolphin 2000). This package uses 
model PSFs to automatically detect, characterise and per- 
form photometry of objects in ACS images. Furthermore it 
incorporates aperture corrections, charge transfer efficiency 
(CTE) corrections (Reiss 2003) and transformation of HST 
flight system magnitudes to the standard Johnson-Cousins 
magnitude system (Sirianni et al 2005). Objects classified 
as being good stellar candidates were retained from the 
dolphot output, with all other objects, characterised as too 
faint, too sharp, elongated or extended, being discarded. A 
further cut was performed to remove all objects containing 
bad or saturated pixels, and stars which had PSF fits with 
X>2.5. 

In the 2004 observations SN 2002ap was detected with 
a significance of ~ 8<r in the F8I4W image, a significance 
of ~ 5cr in the F555W image, but was not detected in the 
F435W and F625W images down to their 5cr detection lim- 
its. Photometry of the SN and of several nearby stars (see 
Figure [2)| is given in Table [4] Also shown are 5a detection 
limits for each of the HST observations. Although objects 

2 and 3 in Figure [2] appear to be a single object they are 
detected by dolphot as two separate PSFs. Proper scaling 
of the F435W image allows us to verify that at least two 
sources exist, where object 2 appears as a much fainter ex- 
tension to the bright source number 3. From the positions 
measured by dolphot the two sources are separated by ~ 

3 ACS/HRC pixels (C'075). The diffraction limited resolu- 
tion of HST results in stellar PSFs with FWHMs of 0'.'05 
in F435W and 0'.'09 in F814W. The much larger magnitude 
difference and the increased blending of objects 2 and 3 as 
observed in the F8I4W image makes it impossible to distin- 
guish the two sources by eye in this image. Since Figure [2] 
is a combination of the F435W and F814W observations, 
objects 2 and 3 appear to be a single source. 

Further late-time images of SN 2002ap were taken at 
the 4.2-m William Herschel Telescope in late August 2006 
(see Table [T] for details) . B and R band observations with 
exposure times matching those of the pre-explosion CFHT 
images were made using the the Auxiliary Port Imaging 
Camera (AUX), which has a pixel scale of O'.'ll/pixel. The 
reduction and subsequent combining of the images was per- 
formed using standard tasks within iraf. There were 8 x 600s 
exposures taken in the B band on 2006 August 25th. Inspec- 
tion showed the combined image quality to be 0'.'75. A set of 
6 x 600s R band exposures were taken on 2006 August 27th, 
and the quality of the final combined image was 0'.'7. 



3 DATA ANALYSIS 

3.1 Alignment of pre- and post-explosion 
observations 

The site of SN 2002ap was found to lie on chip 8 of the 
CHF12K instrument in all the pre-explosion observations. 
All the images at CFHT were dithered by the observers 
to improve rejection of hot pixels and other CCD defects. 
Hence we registered the images in each filter by applying 
simple linear shifts and then combined them using standard 

3 http:/ /purcell. as. arizona.edu/dolphot/ 
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Figure 2. HST image showing SN 2002ap and several neighbour- 
ing stellar objects. This image is the sum of the F435W (Aug 
2004) and F814W (Jul 2004) HRC observations. Although ob- 
jects 2 and 3 appear to be a single source in this image, they are 
detected as two PSFs separated by 3 ACS/HRC pixels (0'.'075) 
by the PSF-fitting photometry package DOLPHOT. See text for 
further details. 



tasks within iraf. The stacked BVR images were then co- 
aligned using the i?-band image as the reference frame. To 
determine the precise location of SN 2002ap on the CFHT 
pre-explosion images, we employed the same method as used 
by, for example, Smartt et al. (2004) and Maund & Smartt 
(2005). 

In order to calculate a transformation to map the 
ACS /HRC (post-explosion) coordinate system to the CFHT 
(pre-explosion) observations, the positions of stars common 
to both sets of images were required. It proved particularly 
difficult to identify common stars due to the extreme differ- 
ences in image resolution (CFHT ~f/.'8; HST ~0'.'08) and 
the pixel scales of the two instruments (CFH12K O'.'2/pixel; 
ACS/HRC O'.'025/pixel). Most of the stars in the HRC field 
of view are in resolved clusters, which appear blended in the 
lower resolution ground based observations. Such stars can- 
not be used to align the images since reliable measurements 
of their positions cannot be made in the CFHT frames. Nev- 
ertheless ten isolated stellar objects common to both the pre- 
and post-explosion observations were identified, and their 
positions measured using the centring algorithms within the 
iraf daophot package. 

A general geometric transformation (x and y shifts, 
scales and rotations) which mapped the ACS/HRC coor- 
dinate system to that of the CFH12K instrument was cal- 
culated using the iraf task GEOMAP and the positions of 
those stars common to both sets of observations. Figure [3] 
shows sections of the pre- and post-explosion images centred 
on the SN coordinates and aligned using the GEOMAP solu- 
tion. The pixel position of SN 2002ap in the post-explosion 
F8I4W frame was measured using the three centring algo- 
rithms of the daophot package (centroid, ofilter, gauss), 
and the mean value transformed to the coordinate system 
of the CFHT images using the iraf task geoxytran. The 
error in the SN position was estimated from the standard 
deviation of these three measurements to be ~ 0'.'003. This 
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Figure 3. Pre- and post-explosion imaging of the site of SN 2002ap. All images are centred on the SN location. The top panel is the 
pre-explosion CFHT R band observation, the middle panel is the post-explosion HST F625W observation (120 seconds) taken in Jan 
2003, and the lower panel is the sum of the post-explosion HST F814W (720 seconds) and F435W (840 seconds) observations taken in 
Jul/Aug 2004. Note that the F625W image is much shallower than either of the constituent observations that make up the bottom panel, 
which explains why it is difficult to identify objects common to both these images. The F625W observation is shown here to convince 
the reader of the location of SN 2002ap. The SN has faded significantly between the 2003 and 2004 observations, but as can be seen from 
Figure [2] it is still visible at the later epoch. 
© 2007 RAS, MNRAS 000, 000-000 
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Table 4. Photometry of SN 2002ap and nearby stars in the HST observations of 2004 (see Figure [2j. Transformation from HST flight 
system magnitudes to Johnson-Cousins magnitude system made using the transformation of Sirianni et al. (2005). Also shown are the 5<r 
detection limits derived for each filter at the location of the SN. The detection limits were transformed to standard magnitudes assuming 
zero colours. 





B 


V 


R 


I 


5cr limit 


26.3 


25.7 


25.4 


25.6 


1 


25.43(0.10) 


23.48(0.04) 


22.16(0.03) 


21.45(0.01) 


2 


24.49(0.06) 


24.60(0.09) 


24.86(0.18) 


24.87(0.17) 


3 


22.72(0.02) 


22.06(0.02) 


21.72(0.02) 


21.45(0.01) 


4 


24.85(0.07) 


25.29(0.13) 


24.89(0.15) 


25.22(0.16) 


5 


25.30(0.09) 


25.56(0.16) 






SN 




25.84(0.19) 




25.05(0.13) 



Figures in brackets are the photometric errors 




Figure 4. The original pre-explosion observation (left), and the resultant image after PSF fitting and subtraction of objects close to the 
SN site (right). 



error proved insignificant when compared to the RMS error 
of the transformation, which was ~ 0'.'072. 



this, two methods were used to determine whether or not 
the progenitor star had been recovered. 



3.2.1 



3.2 The SN site in pre-explosion observations 

The location of the SN site in the CFHT images was inves- 
tigated to discover if a progenitor star had been observed. 
Figure [3^ shows that the SN position lies on the edge of a 
bright source, which is blended in the CHFT frames but, 
as can be seen from the HST imagery (Figure [5}, is formed 
from at least five separate objects. The FWHM of stellar 
sources in the CFHT BR frames is ~ 0'.'8 and in V ~ l'.'2. 
The SN was found to lie some 2" from the nearest of the 
sources that constitute the blend, so any progenitor star in 
the CFHT images should be resolved, provided the observa- 
tions were deep enough for it to be detected. 

Visual inspection showed there to be significant flux 
close to the SN location in the B and R frames, with nothing 
visible in the much shallower V band image. At first sight 
this appears to be a detection of a progenitor star. To test 



PSF fitting photometry using DAOPHOT 



Method 1 utilised the PSF photometry tasks within the iraf 
DAOPHOT package, and followed the techniques for crowded 
field photometry described by Massey & Davis (1992). Aper- 
ture photometry was initially performed on objects across 
the CFHT frames, before an empirical model point spread 
function (PSF) was created from several isolated stars using 
the task PSF. Several different PSF models were created in 
order to test how sensitive the photometry results were to 
the model used. The task ALLSTAR was run on a subset of 
stars in the CFHT frames for each PSF model, and the pho- 
tometry results from each run compared. The results varied 
little with the PSF model used, and could be considered 
identical to within the photometric errors. The positions of 
the SN and of the nearby stars were measured in the HST 
images and transformed to the CFHT coordinate system 
using the transformation discussed in Section[3]T] These po- 
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Table 5. DAOPHOT PSF fitting photometry of supernova site in 
pre-explosion CFHT observations. 



Filter Aposn (mas) Mag 

Recentering during profile fitting 

B 435 25.47(0.10) 

R 368 24.80(0.14) 

No recentering during profile fitting 

B 25.85(0.19) 

R 25.00(0.19) 

Figures in brackets are the photometric errors 



sitions, along with the magnitudes measured by dolphot, 
were fed into the ALLSTAR task, which performed simulta- 
neous PSF fitting on all the objects, including any possible 
progenitor at the SN site. The success of the fitting pro- 
cedure was judged from the fit residuals in the subtracted 
image. Figure U shows an original CFHT frame and a PSF 
subtracted version. 

The ALLSTAR task was implemented in two modes: 1) 
allowing re-centring of objects during PSF fitting using the 
initial input coordinates as starting positions, and 2) no re- 
centring of objects. A summary of the results is presented 
in Table [5] 

In the case where re-centring was permitted the position 
of the object recovered closest to that of the SN was signif- 
icantly shifted from the supernova's precise location; 0'.'435 
in the B band image and 0'.'368 in the R band. The posi- 
tions of this object in the pre-explosion observations were 
measured using several methods and mean values were cal- 
culated for each filter. One measurement was made by fit- 
ting all of the sources simultaneously. Three more measure- 
ments were made by first of all subtracting out all of the 
sources except the possible progenitor, and then measur- 
ing its position using the three centring algorithms of the 
iraf DAOPHOT package. One final measurement was made 
by fitting a PSF to the source in this subtracted image. The 
positional errors were estimated from the standard devia- 
tion of these five values and are shown in Table [6] along 
with the other astrometric errors discussed in Section 13.11 
From the total astrometric errors appropriate for each filter 
we find 4.1ct and 2.5a differences between the position of 
the progenitor site and the object positions measured in the 
B and the R band images respectively. Moreover, the posi- 
tions measured in the pre-explosion B and R band frames 
are themselves separated by ~ (X'200. The astrometric error 
with which to compare this separation is independent of the 
transformation error, and is simply a combination of those 
associated with measuring positions on each of the CFHT 
frames; 0'.'148. Therefore 0'.'200 is equivalent to a separa- 
tion of ~ 1.4a, which is too high to confirm coincidence 
yet too low to rule it out. Equally the displacements of the 
pre-explosion sources from the SN position are of such signif- 
icance that it is impossible to draw any definite conclusions 
on coincidence. 

However it is possible that the re-centred positions of 
these objects are unreliable. It has been found that while 
allowing stars to be re-centred during profile fitting, the po- 



Table 6. Astrometric error associated with the alignment of pre 
and post-explosion observations. 



Error Value (mas) 

SN position 3 

Geometric transformation (RMS) 72 

Pre-explosion position - B 78 

Pre-explosion position - R 126 

Total Error - B 106 

Total Error - R 145 



sitions of some of the faintest stars may be greatly shifted 
towards peaks in the subtraction noise of nearby very bright 
stars (Shearer et al. 1996). Since no objects corresponding 
to those found in the CFHT images were observed in the 
HST frames, and since the shift from the SN position was 
towards the nearby bright stars, it is reasonable to suggest 
that the above might have happened in this case. 

The alternative is therefore to prohibit re-centring, forc- 
ing each PSF fit to be centred at the input coordinates 
throughout the fitting process. Via this method it was possi- 
ble to fit a PSF at the progenitor site in both the B and the 
R band images. The magnitudes of these fits were slightly 
fainter than those from the re-centring case; ~ 0.35 mags 
fainter in B and ~ 0.20 mags fainter in R (Table [5}. How- 
ever, by disabling re-centring we lost any independent verifi- 
cation of the object's position, since we forced the PSF fit at 
the predetermined SN position. It is therefore not possible 
to conclude from PSF fitting techniques that the progenitor 
star of SN 2002ap is detected in the pre-explosion images, al- 
though the fact that the visible flux cannot be explained by 
any objects in the HST post-explosion observations makes 
such a conclusion appealing. 



3.2.2 Comparison with ground based late-time 
observations 

Another method that can be used to detect the presence of 
a progenitor star in archival imagery is to subtract a late- 
time post-explosion image from the pre-explosion version. 
The progenitor star will obviously not exist in the post-SN 
images, and as long as the SN itself has faded enough to 
contribute negligible flux, the subtraction of post-SN im- 
agery from pre-SN frames should reveal a positive source at 
the progenitor position. This of course depends on whether 
the original observations were deep enough for the progen- 
itor star to be detected above the background noise. To 
perform such an image subtraction late-time observations 
of SN 2002ap were taken with the Auxiliary Port Imag- 
ing Camera (AUX) of the 4.2-m William Herschel Telescope 
(WHT) (details in Table [TJ. Observations were made using 
B and R band filters, with identical exposure times and in 
similar seeing conditions to the CFHT observations. Taken 
approximately four and a half years post-explosion, the SN 
was much too faint to have been detected. (Note the mag- 
nitudes and detection limits for the SN in the 2004 HST 
observations - Table [4]) Visual inspection showed the depth 
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of the pre- and post-explosion images to be very similar (see 
Figures [5^, and [5p) . 

After re-binning the WHT frames to match the pixel 
scale of the CFH12K camera (CFH12K pixel scale 072; AUX 
pixel scale O'.'l) and accurate image alignment, the image 
subtraction package ISIS 2.2 (Alard & Lupton 1998; Alard 
2000) was used to subtract the post-explosion images from 
their pre-explosion counterparts. Subtraction of the B-band 
observations resulted in a near perfectly fiat subtracted im- 
age at the position of the SN (see Figure [5}i), implying that 
the same source is visible pre- and post-explosion and there- 
fore cannot be the SN progenitor. However, in the case of the 
B-band subtraction an extended negative feature, albeit of 
quite low signal-to-noise, is visible overlapping the SN po- 
sition (see Figure [BJ;). Given the sense of the subtraction 
(pre minus post) a negative feature indicates an increase in 
detected flux in the post-explosion B-band image, an obser- 
vation that is inconsistent with the detection of a progenitor 
in the pre-explosion frame. Also visible in the B-band sub- 
tracted image is a bright positive residual coincident with 
the centre of the adjacent blended source some 2" from the 
SN position, implying that something in this blend appeared 
brighter in the pre-explosion image. We can instantly rule 
this out as a progenitor detection based solely on the astrom- 
etry. One might argue that the progenitor is detected in the 
pre-explosion frames and that the flux in the WHT post- 
explosion B and B-band observations is from the SN, which 
has rebrightened perhaps due to interaction with circum- 
stellar material or a light echo from interstellar dust. This 
scenario could potentially explain the differences between 
the B-band images. On the other hand it is highly unlikely. 
Comparing the B and R magnitudes of the pre-explosion 
source (Table [5} with the 5a detection limits of the SN in 
August 2004 (Table [U we see that by this epoch the SN has 
already become significantly fainter than the pre-explosion 
object. In order to reproduce the subtracted images, the SN 
would have to become significantly brighter between August 
2004 and August 2006 when the WHT images were taken. 
Crucially, its R magnitude at this later epoch would have to 
exactly match that of the pre-explosion source, with its B 
magnitude becoming only slightly brighter. Such coinciden- 
tal matching of the SN magnitudes to those of the pre-SN 
source is virtually impossible. 

A simpler explanation is that the source close to the SN 
position in the pre- and post-explosion images is the same 
object, with the residuals in the subtracted image arising 
from differences between the CFH12K and the AUX filter 
functions. Since this object is still visible post-explosion ob- 
viously it cannot be the progenitor. So what is it? 

Nothing is visible in the 2004 HST observations at the 
position of this source, which lies between the SN and the 
closest of the nearby stars in Figure [2] This is in spite of the 
fact that the detection limits of the HST images are signifi- 
cantly deeper than the measured magnitudes of this object 
in the CFHT frames. If, however, we assume that the object 
is a region of extended emission, we can explain this non- 
detection by the apparent low sensitivity of the ACS/HRC 
to such sources when compared with the CFHT instrument. 
Inspection of Figure [3] shows several areas where there is sig- 
nificant flux in the CFHT images which cannot be accounted 
for by stars visible in the high resolution HST frames. Some 
traces of this flux are just visible in the HST F814W obser- 



vation from 2004. Due to the very small pixel scale of the 
ACS/HRC (O'.'025/pixel) the light from a region of extended 
emission is spread over many pixels and is therefore easily 
dominated by detector read noise if exposure times are not 
sufficiently long. In contrast the CFH12K pixels, which cover 
an area on the sky sixty four times that of the HRC pixels, 
are much more sensitive for imaging such sources. Summing 
the HST frames from all filters, and rebinning to match the 
pixel scale of the CFH12K camera, gave some indication of 
a diffuse region of emission close to the SN position and 
extending towards the bright stars nearby, but attempts to 
estimate its magnitude proved futile due to the very low 
signal-to-noise ratio. 

If the flux from this diffuse source is dominated by emis- 
sion lines (for example a Hn region) this fact coupled with 
the differences between the CFHT and the WHT filter func- 
tions could explain the negative residuals in the subtracted 
B-band image (see Figure [5} . Unfortunately the Ha pre- 
explosion observation was too shallow (900s compared to the 
3600s B-band exposure) to produce a significant detection 
at the SN site, so it cannot be used to confirm the presence 
of a H II region. 

We conclude that the object seen in the pre-explosion 
observations (in Fig. [3} close to the SN position is a diffuse 
source, possibly a Hn region, and not the SN progenitor. 

3.3 Detection limits of the pre-explosion 
observations 

Since the progenitor star was not detected, limiting magni- 
tudes were derived for the CFHT B and R band observations 
using several methods. 

The standard deviation of the counts per pixel, in an 
area of 5x5 pixels centred on the SN location, was used as a 
measure of the noise per pixel at the SN site in each image. 
By rearranging the equation for Signal-to-Noise ratio 

S/N = Fstar (1) 

where Fstar is the flux of a star and a is the background 
noise level (both in electron counts), one can calculate what 
Fstar must be to produce a detection with a signal-to-noise 
ratio of a given value. In this case we decided to calculate 
the Fstar that would give a 5a detection within a 4 pixel (« 
FWHM of the PSF) radius aperture. An aperture correction 
was applied to the resultant magnitude to give the final 5a 
limiting magnitude. The value of a in Equation [1] is the 
standard deviation in electron counts for the whole 4 pixel 
aperture, a ap , and is calculated by multiplying the measured 
standard deviation per pixel, a P i X , by the root of the number 
of pixels in the aperture, \J~N. 

Another method involved measuring a ap directly. A 6x6 
grid of 4 pixel radius apertures was used to measure the flux 
of an area of blank sky adjacent to the SN location. For 
each aperture the sky background was measured in a con- 
centric annulus and subtracted from the total counts in the 
aperture. a ap was then calculated as the standard deviation 
of the remaining counts in the 36 apertures. This method 
is preferable to using the total counts in each aperture to 
calculate the standard deviation since total counts can be 
significantly affected by changes in the level of the smooth 
background. An aperture at one side of the grid may overlie 
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Figure 5. a) CFHT pre-explosion B observation, b) WHT post-explosion B observation, c) pre minus post subtracted S-band image 
and d) pre minus post subtracted i?-band image. All the images are shown with inverted intensities (positive sources appear black, 
negative sources appear white). The ij-band subtraction (d) appears to be perfectly flat at the position of the SN, which is marked in 
all frames by cross hairs. In the f?-band subtracted image (c) a faint and extended negative (white) feature is visible overlapping the SN 
location. Adjacent to this, coincident with the centre of the blend of nearby stars, is a bright positive (black) source. It is believed that 
both of these subtraction residuals are due to differences in the transmission functions of the CFHT and WHT filters (see Figure [6jl, 
where the negative feature overlapping the SN position is a diffuse source of emission line flux (possibly a Hn region) that is more 
efficiently tranmitted by the WHT B filter, and the positive source is due to extra continuum flux from the nearby bright stars that is 
transmitted by the CFHT filter. 



a region of much higher background flux than an aperture at 
the other side, producing large count differences and hence 
a value of a ap that is unrealistically high. Subtraction of the 
local background at each location produces aperture counts 
that are more readily comparable, as they are measures of 
the deviation from an average level. Since we wished to mea- 
sure the background fluctuations and noise on a scale com- 
parable to the size of the PSF in the images and not large- 
scale variation in the background it was considered best to 
use the background subtracted counts. 

One final estimate for the 5<r detection limits was made 
using the photometry of stars within a 100" x 100" section 
of each image, roughly centred on the SN position. The 5a 
detection limit was determined as the average magnitude of 
all the stars with a photometric error of ~ 0.2 mags. 



The detection limits determined from the above meth- 
ods were averaged to produce 5<r limiting magnitudes of B = 
26.0 ±0.2 mags and R = 24.9 ±0.2 mags. A magnitude limit 
for the much shallower and less restrictive V band image 
was not derived. 



4 DISTANCE, EXTINCTION AND 
METALLICITY 

Hendry et al. (2005) used three different methods to esti- 
mate the distance to M 74; standard candle method (SCM), 
brightest supergiants, and a kinematic distance estimate. 
From the mean of these three methods they estimate the dis- 
tance to be 9.3±1.8 Mpc. A light echo was discovered around 
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another SN which exploded in this galaxy (SN 2003gd) and 
this study suggests a smaller distance of 7.2 Mpc would 
be more consistent w ith a model for the light echo flux 
jVan Dvk et al. 20061 ). A further determination of distance 
from the expan ding photosphere method gives a distance of 
6.7 ±4.5 Mpc (|Vinko et al. 20 02). However given the large 
uncertainty of this result and the fact that the light echo 
method does not produce a definite independent estimate, 
we will adopt the distance of 9.3 ± 1.8 Mpc, noting that 
if the distance is lower the luminosity limit of the progeni- 
tor would be closer to the lower limit of the uncertainty we 
derive. 

The total extinction toward SN200 2ap has been shown 
by all previous s tudies to be low (Smartt ct al. 20021 ; 
IVinko et al. 20021 ; iFolev et al. 20031 ; iMazzali et al. 20021 ) 
and we adopt the value of E(B — V)= 0.09 ± 0.01 measured 
from a high resolution spectrum of the SN by Takada-Hidai 
et al. (2002). Assuming the reddening laws of Cardelli et al. 
(1989), with R v = 3.1, gives A B = 0.37 ± 0.04 and A R = 
0.21 ± 0.02. 

The metallicity gradient of M74 has been studied from 
the nebular emission lines of HII regions by Van Zee et al. 
(1998) and we have estimated the oxygen abundance from 
the latest emission line strengths combined with the latest 
empirical calibration from Bresolin et al. (2004). The abun- 
dance gradient in M7 4 suggests that at the galactocentric 
distance of SN2002ap |Smartt et al. 20*021 ) the metallicity is 
12 + log [O/H] = 8.4 ± 0.1. This latest calibration tends 
to give lower abundances for significantly metal rich regions 
than have previously been derived, but does not significantly 
change that originally derived in Smartt et al. (2002). An- 
other recent study of the metallicity gradient in M74 is that 
of Pilyugin et al. (2004). At the position of SN 2002ap this 
would suggest a metallicity of 8.3±0.1 dex. Hence it appears 
that the environment of SN2002ap is mildly metal deficient, 
somewhat similar to the massive stars in the Large Magel- 
lanic Ck2ud i _wl2ich_liave an oxygen abundance of 8.35 ± 0.1 
dex l|Hunter et al. 2006T ). However we note that Modjaz et 
al. (2007) determine a slightly higher value of 8.6 dex using 
a different calibration of the oxygen line fluxes. 



5 THE PROGENITOR OF SN 2002AP 

From the detection limits of Section 13.31 and the values of 
distance and extinction from Section [4] we derive 5a ab- 
solute magnitude limits for the progenitor of SN 2002ap of 
M B > -4.2 ±0.5 and Mr 5? -5.1 ±0.5. The uncertainties are 
the errors in the sensitivity limits (±0.2), distance modulus 
(±0.43) and extinction (A s ±0.04; Ar±0.02) combined in 
quadrature. Although no progenitor star is detected, these 
limits constitute the tightest constraints yet placed on the 
absolute magnitude of a Type Ic SN progenitor, and as such 
can be used to deduce the likely properties of the star. The 
properties of greatest interest are the type of star that ex- 
ploded and its mass, and several methods were employed to 
estimate these from the data. 

5.1 Single star evolution models and luminosity 
limits 

Since we had no knowledge of the object's colour, we calcu- 
lated 5a luminosity limits from Mb and Mr for a range of 
supergiant spectral types using the colours and bolometric 
corrections of Drilling & Landolt (2000). This method is the 
same as that used by, for example, Maund & Smartt (2005). 
The luminosity limits were plotted on Hertzsprung-Russell 
(H-R) diagrams (Figure [7}, along with stellar evolutionary 
tracks for single stars of initial mass between 7 and 120Mq. 
Four sets of stellar models were used, all of approximately 
LMC metallicity (metal/mass fraction Z=0.008); two sets 
from the Geneva group (Schaerer et al. 1993, Meynet et al. 
1994) and two sets created using the Cambridge STARS3 
code (Eldridge et al. 2006). Of the two sets of models from 
each code, one incorporated standard mass-loss rates and 
the other arbitrarily doubled mass-loss rates during the pre- 
WR and WNL phases of evolution. Models from both groups 
were used in order to gauge any model dependence in our 
interpretation. Any supergiant star belonging to the solid 
shaded regions above the luminosity limits in these H-R di- 
agrams would have been detected in one or more of the 
CFHT images, and can therefore be ruled out as a possible 
progenitor. In this way all evolved hydrogen-rich stars with 

4 http:/ /www. ast.cam.ac.uk/~stars 
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initial masses greater than 8Mq can be ruled out as possible 
progenitors. Such a conclusion is entirely consistent with the 
spectral characteristics (lack of H-lines) of SN 2002ap. 

At effective temperatures (T e g) higher than those of the 
hottest O-stars are the evolved stages of stars with initial 
masses greater than -~25-30Mq. At this late stage of their 
evolution such stars have lost all or most of their H-rich en- 
velope due to mass loss processes, such as a strong stellar 
wind, explosive outbursts or interaction with a binary com- 
panion. The resulting objects (Wolf-Rayet (W-R) stars) are 
the exposed helium cores of the original massive stars. W-R 
stars are divided into two main subtypes; WN stars, in which 
the emission spectra are dominated by lines of nitrogen and 
helium, and WC stars, where carbon and oxygen lines dom- 
inate. The most oxygen rich WC stars are often further clas- 
sified as WO. These sub-types are best explained as an evo- 
lutionary sequence of the form WN — > WC — > WO (Conti 
1982). The enhanced helium and nitrogen abundances ob- 
served in the spectra of WN stars correspond well to the 
equilibrium products of H-burning via the CNO cycle, which 
appear at the stellar surface as a result of mass loss. Sub- 
sequent He-burning and further mass loss results in the en- 
hancement of carbon and oxygen at the stellar surface along 
with the depletion of helium and nitrogen. The star has be- 
come a WC star. Still more advanced stages of helium burn- 
ing result in an overabundance of oxygen at the expense of 
carbon, producing a WO star. 

SN theory predicts that the progenitors of Type Ib/c 
SNe are W-R stars, since Ib(Ic) spectra show no signs of 
H(H or He), which is most easily explained by the progeni- 
tor's lack of a H-rich(He-rich) envelope. A the oretical Mass- 
Luminosity relationshi p exists for W-R stars l|Maeder 19831 ; 
ISmith fc Maeder 19891 ; ICrowther 20071 ). Crowther (2007) 
states this relationship for a H-free W-R star as 

log ± = 3.032 + 2.695 log **- - 0.461 (log |-) " . (2) 

Using this formula one can calculate the mass of a W-R 
star if its luminosity is known. Note this calculation yields 
the mass of the W-R star, not the initial mass of the star 
from which it evolved. In this case we can calculate an up- 
per limit for the W-R star's mass from the sensitivity limits 
of the CFHT observations. However to convert our Mb and 
Mr limits to luminosity limits requires appropriate bolo- 
metric and colour corrections. Since it is not possible to 
constrain the colours of the progenitor star from magni- 
tude limits we instead consider the broadband colours of 
W-R stars in the LMC (Massey 2002), the vast majority 
of which have colours in the ranges —0.4 ^ B ~ V ^ 0.6 
and —0.1 ^ V — R ^ 0.4. Bo lometric corrections range from 
-2.7 to -6.0 (|Crowther 20071 ). and Smith & Maeder (1989) 
find a BC of —4.5 ± 0.2 as appropriate for a wide range of 
galactic W-R stars. The bolometric correction is, however, 
related to the narrowband v magnitude of the star, as pho- 
tometry of W-R stars is usually performed with narrowband 
filters in order to sample the continuum flux between strong 
emission features. Differences between narrow (denoted by 
lowercase letters) and broadband photometry of W-R stars 
can be considerable due to the extra flux contributed by 
emission lines in the broadband observations. For example, 
b — B — 0.55 and v — V = 0.75 have been found for the 
strongest lined WC stars (Smartt et al. 2002). Here we use 



our broadband sensitivity limits and note that the actual 
bolometric magnitude limits will be equal to or fainter than 
what is calculated. For example if 6 — B — 0.55 is applicable 
for the progenitor, then the limiting Mboi should be fainter 
by 0.55 mag. 

Approximate values of Mboi were calculated from the 
magnitude limits Mb and Mr using a BC of -4.5 and the 
median of the colours of the LMC W-Rs appropriate in each 
case. 

M B -> Mboi > -8.80 ± 0.60 
M R -> M bo i > -9.45 ± 0.56 

The errors are a combination of the estimated uncer- 
tainties in the absolute magnitudes (±0.5), bolometric cor- 
rection (±0.2), and colours (B-U±0.26; V-R±0.17). The 
luminosity limit derived from the most constraining value of 
Mboi is log(L/L ) < 5.40±0.24. Using equation [2] this limit 
corresponds to a final W-R progenitor mass of < 12±|M Q , 
which is entirely consistent with that of ~ 5Mq found by 
Mazzali et al. (2002) from modelling of the SN explosion. 

It may be possible to infer a likely range of initial masses 
through comparison with stellar models as we have done for 
supergiant stars. For completeness our W-R luminosity limit 
is plotted in Figure [7] as a straight line across each plot at 
high temperatures. This is a somewhat crude approximation 
but since the error is conservative, the upper limit allows for 
equally conservative constraints to be placed on the star's 
properties. Considering the models with standard mass loss 
rates first, the Geneva and the STARS codes predict W-R 
stars to form only from stars initially more massive than 
about 34M0. These independent stellar evolution codes ap- 
pear to be in good agreement, although there are noticeable 
differences during W-R evolution phases. These differences 
are manifested in the initial mass versus final mass functions 
produced by each code, the final luminosity of each model 
W-R star being representative of final mass. The STARS 
models produce higher final W-R masses with increasing 
initial mass. The Geneva tracks show a similar behaviour, 
but reach a peak final mass at an initial mass of 6OM0 . The 
Geneva models of initial mass higher than 60M© produce W- 
R stars of substantially lower final mass than similar STARS 
models. These differences are due to the lower mass loss rates 
prescribed during W-R evolution in the STARS code, which 
uses the W-R mass loss rates of Nugis & Lamers (2000) and 
includes scaling of W-R mass loss rates with initial metal- 
licity (Eldridge & Vink 2006) . Despite the discrepancies the 
overall conclusions drawn from both sets of standard mass 
loss models are the same. Ultimately at no point during W-R 
evolution do any of the models fall below the detection lim- 
its. This implies that no singly evolving star at metallicity 
Z=0.008 could have been the progenitor of SN 2002ap. 

If single-star evolution is to be used to explain the pro- 
genitor of SN 2002ap, then we must invoke enhanced mass- 
loss in the models. As mentioned previously, mass loss can be 
a result of several phenomena such as stellar winds, eruptive 
outbursts and binary interaction. The latter we consider in 
Section T5. 31 For now we consider only singly evolving stars. 
Eruptive events such as LBV bursts are poorly understood 
but may account for substantial amounts of mass lost by 
very massive stars during their lifetimes (Smith & Owocki 
2006). Wind-driven mass- loss is slightly better understood. 
Three main factors affect this process: luminosity, surface 
temperature and composition. Higher luminosities in turn 
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Figure 7. H-R diagrams showing the 5<r luminosity limits for the progenitor of SN 2002ap and stellar evolutionary tracks for singly 
evolving stars with initial masses of 7-120Mq. Top left and top right show the Geneva and Cambridge STARS stellar models respectively, 
both for metallicity Z=0.008 (approximately LMC) and incorporating standard mass loss rates. Bottom left and bottom right show the 
Geneva and STARS models at the same metallicity but with double the standard mass loss rates during the pre-Wolf-Rayet and WNL 
phases of evolution. Any star within the shaded area on each plot would have been detectable in one or both of the BR pre-explosion 
observations. Only those models falling outside of this shaded area can be considered as viable progenitors for SN 2002ap. 



drive stronger stellar winds and increased metallicity also 
leads to higher mass loss rates (e.g. Abbott 1982). This re- 
sults in the minimum initial mass for a W-R star decreas- 
ing with increasing metallicity. The Geneva group provide 
Z=0.008 models with double the standard mass loss rates 
during the pre- W-R and WNL phases of evolution (i.e. while 
each model retains some hydrogen), and we produced simi- 
lar STARS models for comparison. The increased mass loss 
rates may be considered as a correction for an underesti- 
mation of the rates at Z=0.008, or of the initial metallicity 
of the progenitor of SN 2002ap. A factor of two increase 
in mass-loss rates is believed to represent a realistic upper 
limit, given the intrinsic uncertainty of the rates and the 
uncertainty in our measurement of metallicity. Comparing 
the luminosity limit to these tracks provides robust limits on 
what types of single stars, evolving in the canonical fashion, 
might reasonably have been the SN progenitor. The Geneva 
tracks with x 2 standard rates suggest the progenitor would 
have to have been initially more massive than AQMq. Tak- 
ing into account the upper error on the luminosity limit, the 
STARS x2 models suggest that the progenitor could have 
evolved from a single star with initial mass in the range 
30-60Mq, or a star of 12OM0. The differences between the 
Geneva x2 and STARS x2 models are again due to the 
different W-R mass-loss rates employed by each code. 



It is interesting to compare our analysis here with the 
work of Maund, Smartt & Schweizer (2005). In an attempt 
to place limits on the progenitor of another Type Ic su- 
pernova, SN 2004gt, they calculated colour and bolometric 
corrections for W-R stars by performing synthetic photom- 
etry on model W-R spectra from Grafener et al. (2002) (see 
Maund & Smartt 2005 for details). In this manner they de- 
rived a typical bolometric correction of ~ -2.7. Applying a 
similar BC to the limiting magnitudes presented in this pa- 
per would push our W-R luminosity limit down by 0.7 dex 
to a median value of log(L / '£©) = 4.7. The upper W-R mass 
limit inferred in this case would be around 5Mq , much more 
restrictive than the 9Mq limit found by Maund, Smartt & 
Schweizer (2005) for the progenitor of SN 2004gt. This re- 
vised luminosity limit would rule out all single star progeni- 
tors, even from those models with double the standard mass 
loss rates. Our choice of a BC of -4.5 has been explained 
above and although the limits produced using this value are 
conservative in comparison, we believe they offer the most 
realistic interpretation given the model uncertainties. 

In Figure [8] we compare our Mb and Mr limits with 
broadband photometry of W-R stars in the LMC (Massey 
2002), a comparison that is independent of BC. Considering 
the more restrictive Mb limit, around 79 per cent of WNL 
stars can be ruled out as possible progenitors along with 
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23 per cent of WNE stars and 75 per cent of WC and WO 
stars. Given that SN 2002ap was a Type Ic SN, and therefore 
deficient of both hydrogen and helium, one would expect the 
W-R progenitor to have been an evolved WC or WO star. 
Assuming this to be true we can say that the progenitor was 
amongst the faintest ~ 25 per cent of WC/WO stars which 
are known to exist. However our sensitivity limits cannot 
directly rule out around 21 per cent of the WNL stars and 
77 per cent of the WNE stars as possible progenitors. WNL 
stars, which mark the transition of normal stars to W-R 
stars, have not yet lost all of their hydrogen envelope and can 
almost certainly be ruled out as progenitors based on the SN 
classification. Since such a large fraction of observed WNE 
stars fall below our detection limits, and with uncertainty 
as to how much helium must be present to result in a Type 
lb instead of a Type Ic SN, we cannot definitely rule out a 
WNE star progenitor. 

5.2 Combining observational constraints 

Besides trying to constrain the properties of the progenitor 
of SN 2002ap through direct observation, there are various 
indirect methods which can be used. One such approach is 
to observe the SN itself. Through modelling the explosions 
of carbon-oxygen (C+O) stars Mazzali et al. (2002) created 
best fits to the lightcurve and spectra of this SN, and thereby 
estimated that around 2.5Mq of material was ejected dur- 
ing the explosion. Including the mass of the compact stellar 
remnant suggests a C+O progenitor star mass of ~5Mq, 
which would form in a He core of ~7Mq. To produce a 
helium core of this mass requires a star of initially greater 
than ~15-20Mq. (This initial mass estimate is lower than 
the 20-25M Q suggested by Mazzali et al. (2002) due to the 
inclusion of convective overshooting in the STARS stellar 
models.) 

X-ray and radio observations of the SN can provide 
some constraints. The interaction of the fast moving super- 
nova ejecta with the circumstellar material (CSM) leads to 
emission at X-ray and radio frequencies. Interaction with a 
higher density of CSM results in stronger emission. However, 
the radio and X-ray fluxes detected for SN 2002ap soon after 
explosion were low (Berger et al. 2002; Sutaria et al. 2003; 
Soria et al. 2004; Bjornsson & Fransson 2004) suggesting a 
low density for the material immediately surrounding the 
star. If we assume that the final mass of the progenitor was 
around 5Mq, and that it was initially more massive than 
15-2OM0, then the star must have lost 1O-15M0 during its 
lifetime. This material must be sufficiently dispersed prior 
to the SN explosion so as not to result in strong radio and 
X-ray emission. Only two epochs of X-ray observations of 
the SN are available, the first taken less than five days post- 
explosion and the second about one year later (Soria et al. 
2004). Radio observations were taken on a total of 16 epochs, 
15 of which spanned the first ten weeks after discovery of the 
SN (Berger et al. 2002; Sutaria et al. 2003) and one final ob- 
servation approximately 1.7 years post-explosion in which 
the SN was not detected (Soderberg et al. 2006). It is rea- 
sonable to assume that no interaction occurred between the 
SN ejecta and any dense CSM during the 18-month gap in 
the radio monitoring, since such interaction would still have 
been detectable in the late-time observation (see Fig. 2 of 
Montes et al. (2000) - the timescale for variation in the radio 



light curve of SN 1979C is of the order of 4 years). We can 
say that for at least the first 1.7 years after the explosion of 
SN 2002ap the ejecta encountered no region of dense CSM. 
By multiplying this time period by the velocity of the out- 
ermost layers of the SN ejecta we can calculate a lower limit 
for the radius of such CSM. Mazzali et al. (2002) derive 
a photospheric velocity of ~30,000 kms -1 from a spectra 
taken just 2 days post-explosion (Meikle et al. 2002) . How- 
ever, Mazzali et al. (2002) also suggest that the large degree 
of line blending in this spectrum requires sufficient material 
at velocities > 30,000 kms -1 , and subsequently employ an 
ejecta velocity distribution with a cut-off at 65,000 kms -1 in 
their models. Furthermore, through modelling of X-ray and 
radio observations, Bjornsson & Fransson (2004) derive an 
ejecta velocity of ~ 70, 000 kms -1 . They also point out that 
the photospheric velocity only provides a lower limit to the 
velocity of line-emitting regions (see Section 4 of Bjornsson 
& Fransson (2004)). From this range of velocities of 30,000 
to 70,000 kms -1 , we calculate the minimum radius for a re- 
gion of dense CSM to be in the range ~ 1.6 to 3.8 x 10 12 km 
(~ 11,000 to 25,000 astronomical units (au)). 

One might also attempt to constrain the properties of a 
SN progenitor by investigating the properties of other stars 
close to its position. If such stars are found to be coeval with 
each other and are close enough to the SN site (within a few 
parsecs) to be considered coeval with the progenitor, then 
the most evolved of these stars are likely to have had initial 
masses similar to the SN progenitor. However, in this case 
SN 2002ap and its closest stellar neighbours are separated 
from each other by distances of around 50 to 100 pc, as 
can be seen from the HST observations (see Figure [2j . We 
cannot reliably assume that that these stars are coeval with 
each other or with the SN progenitor, and therefore cannot 
use their properties to place any strong constraints on the 
progenitor star. 

Nevertheless photometry was performed on these stars 
as detailed in Section [2~2l and the results are presented in Ta- 
ble[4] Absolute BVRI magnitudes of each star were plotted 
against wavelength and compared with spectral energy dis- 
tributions (SEDs) of supergiant stars of spectral type 09 to 
M5 (Drilling & Landolt 2000) to determine an approximate 
T e ff and appropriate bolometric correction for each star. The 
details of how objects 2 and 3 are spatially resolved, at least 
in the F435W image, have already been discussed fully in 
Section [2.21 However, while fitting SEDs to each of the ob- 
jects, it was found that no single SED could be fitted to 
object 3. This may have been due to increased blending of 
objects 2 and 3 in observations at longer wavelengths, which 
in turn would have lead to unreliable photometry, or it could 
indicate that object 3 is in fact the convolution of several 
sources. Assuming the latter scenario and the simplest so- 
lution therein, i.e. that object 3 is a blend of two stars, we 
created a program to iterate over all possible combinations 
of pairs of supergiant SEDs to find the best fit to the data. 
In this way we tentatively estimated the spectral types and 
magnitudes of stars 3a and 3b. All the stars were plotted on 
an H-R diagram and compared to stellar models to estimate 
their masses, which range from 15 to 40M© (see Figure [9}. 
(Note that this range of masses is not affected if we ignore 
3a and 3b.) The lifetime of a 15M© star is of the order of 
15 Myrs, and for a 40M© star around 5 Myr. If all the stars 
formed at about the same time one would expect to find 
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Figure 8. Broadband photometry of W-R stars in the LMC (Massey 2002) compared to the absolute B (bottom) and R (top) magnitude 
limits. Note that the evolutionary sequence for W-R stars is WNL — » WNE — > WC — > WO. 
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Figure 9. H-R diagram of stars close to the site of SN 2002ap and STARS stellar evolutionary models of massive stars at metallicity 
Z=0.008. Since star 3 is not well fitted by a single SED, this source is better reproduced as the convolution of a double star (3a and 3b). 
See text for details. 



that the most evolved stars are also the most massive. How- 
ever the two most evolved stars in Figure [5] appear to be at 
the lower end of this mass range, at around 15-2OM0. This 
could suggest that the higher mass stars are unresolved mul- 
tiple objects rather than single stars, but it is more likely a 
confirmation that these stars are indeed not coeval. 

In Section 15.11 we found that any single star progeni- 
tor would have to be initially more massive than 3O-4OM0 
and that mass loss rates would have to be double the stan- 
dard for the measured metallicity. The final masses of such 
stars given by the two stellar evolution codes are between 10 
and 12Mq. This is significantly larger than the 5Mq esti- 
mate of Mazzali et al. (2002) , but consistent with the mass 
we derived in Section 15.11 using the W-R mass-luminosity 
relationship (not surprising considering both estimates are 
based on our photometric limits). Enhanced mass- loss rates 
due to continuum driven outbursts (Smith & Owocki 2006) 
might play an important role in the evolution of any pro- 
genitor star that is initially more massive than AQMq and 
could help to produce lower mass W-R stars. 

Another possibility is a progenitor produced by the ho- 
mogeneous evolution (Yoon & Langer 2005) of a 12-20M© 
star. If the rotation of a main sequence star is fast enough, 
the rotationally induced mixing occurs on a much shorter 
timescale than nuclear burning. Fresh hydrogen is continu- 
ally mixed into the core until the entire hydrogen content 
of the star is converted into helium. In this way the main 
sequence star transitions smoothly to become a helium star 



of the same mass. While Yoon & Langer (2005) do not con- 
sider LMC metallicity objects if we extrapolate their results 
it is reasonable to expect that the most rapidly rotating 
stars with initial masses in the range 12-20 Mq will lead to 
a Type Ic progenitor in the mass range implied by our lumi- 
nosity limit and the results of Mazzali et al. (2002). Detailed 
models of rapidly rotating massive stars would be required 
to test this hypothesis. 



5.3 Binary evolution 

Binary evolution provides a possible alternative to the sin- 
gle star models discussed previously. Such a scenario is sug- 
gested by Mazzali et al. (2002) for the progenitor of this 
SN, where a star of initially 15-20Mq loses much of its mass 
through interaction with a binary companion, to become a 
C+O star of around 5Mq. 

Using our luminosity limits, it is possible to constrain 
the binary systems that might have produced the progenitor 
of SN 2002ap. We can place limits on the luminosity and the 
mass of the binary companion star in the same way as we 
did for the progenitor. In the case of the companion, how- 
ever, we can use both the pre- and the post-explosion obser- 
vations for this purpose. A companion star was not seen in 
either set of images, but the post-explosion observations pro- 
vided slightly more restrictive luminosity limits, and these 
are plotted on an H-R diagram in Figure [lOl Whereas with 
the SN progenitor we could assume that the star must have 
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Figure 10. H-R diagram showing the 5<r luminosity limits for a 
binary companion to the progenitor of SN 2002ap. These limits 
were derived from the post explosion HST F425W and F814W 
observations. Also plotted are STARS stellar evolutionary tracks 
for singly evolving stars with initial masses of 7-120Mq, at metal- 
licity Z=0.008 and with double the standard mass loss rates. Any 
star within the shaded area would have been detectable in one or 
both of the aforementioned HST observations. If the progenitor 
star did have a binary companion it must have been either a main 
sequence star of < 20Mq or a dormant neutron star/black hole. 



been in an evolved state just prior to exploding, we cannot 
make the same assumption about any possible binary com- 
panion. It is conceivable that such a star could be at any 
stage in its evolution. Comparing our luminosity limits with 
stellar models we can at least set discrete mass limits for the 
companion star at various evolutionary phases. 

Considering the case where the binary companion is ini- 
tially less massive than the SN progenitor, and therefore less 
evolved at the time of the explosion, we see from Figure [10] 
that it must have been a main sequence star of < 20M©. 
Stars more massive than this would have been visible at all 
stages of their evolution. Note that this constraint is on the 
mass and evolutionary phase of the companion at the time 
of the SN; that is after it has accreted any material from 
the progenitor. In the case of a star initially more massive 
than the progenitor, it would have exploded as a SN prior 
to SN 2002ap leaving a neutron star or black hole as the 
binary companion^. This would of course require the binary 
to remain bound after the first SN. Brandt & Podsiadlowski 
(1995) predict that just 27 per cent of binaries consisting of 
two high mass stars will remain bound after the initial SN 
explosion, and of these bound systems 26 per cent will im- 
mediately experience strong dynamical mass transfer lead- 
ing to their merger. Therefore around 20 per cent of high 
mass binaries might be expected to result in a bound sys- 
tem with a neutron star or black hole component, with only 
a fraction configured in such a way as to result in the binary 
interaction required to produce the progenitor of SN 2002ap. 
Such systems are therefore much rarer than those with main 
sequence companions. 

Using the conclusions of Mazzali et al. (2002) as best 
estimates of the initial and final masses of the progenitor, 

5 Nomoto et al. (1994) suggest that an initially more massive 
companion star may also become a white dwarf in some cases. 



we can say that it lost 10-15Mq of material before explod- 
ing. A substantial fraction of this material must have been 
ejected from the binary system rather than accreting onto 
a main sequence companion, otherwise it could have gained 
enough mass to become detectable in the observations (e.g. 
Maund et al. 2004). We can therefore conclude that the bi- 
nary mass transfer process must have been non-conservative 
in cases where the companion was a main sequence star. 
Furthermore, this non-conservative mass transfer must have 
occurred early enough in the progenitor star's lifetime to 
allow the ejected material sufficient time to disperse to a 
large radial distance. This condition must be met to re- 
main consistent with the X-ray and radio observations of 
the SN discussed in Section 15.21 For a black hole compan- 
ion there is no such requirement that mass transfer must be 
non- conservative since a more massive black hole is no more 
visible than a lower mass counterpart. The rate at which a 
black hole can accrete material is however limited (Narayan 
& Quataert 2005) and in many cases the process will there- 
fore be non- conservative. 

Kippenhahn & Weigert (1967) classified binary systems 
which undergo mass transfer into three types, Cases A, B 
and C. We shall first discuss each case assuming a main 
sequence companion of 20Mq. Case A mass transfer is 
where the donor star begins to transfer mass while still on 
the main sequence. This occurs in very close binary systems 
with periods of the order of a day. However, the expansion 
of the star on the main sequence, which causes it to fill its 
Roche lobe and dictates the rate of mass transfer, occurs on 
the nuclear timescale. The companion star has no problem 
adjusting to accept the material at such a low rate, with the 
result that Case A mass transfer is in general a conserva- 
tive process. Our previous argument that conservative mass 
transfer would result in a main sequence companion becom- 
ing visible in our observations therefore suggests that Case A 
mass transfer is not a viable option. Furthermore, while the 
donor star is the more massive, conservative mass transfer 
will result in the shrinking of the binary orbit. Since Case A 
systems are already initially very close, this will often lead 
to the merger of the two objects. Where the merger is of two 
main sequence stars the result will be a single star of mass 
roughly equal to the sum of its components. This object will 
then continue its evolution as a single star, for which we have 
already discussed progenitor scenarios in Section \5. II 

Case C mass transfer is where the donor begins to trans- 
fer mass at the end of core helium burning when it ex- 
pands rapidly, ascending the giant branch for the second 
time. Case C therefore occurs in wider binaries with or- 
bital periods of around 100 days or more, ensuring that the 
donor must expand to become a supergiant before filling 
its Roche lobe. Mass loss will occur at the thermal or pos- 
sibly even the dynamical timescale of the donor star as it 
rapidly expands. The rate at which a companion can ac- 
crete this material will be dictated by its thermal timescale. 
A main sequence companion of lower mass will have a long 
thermal timescale compared to the donor star, and there- 
fore will only accrete material at a fraction of the rate at 
which it is being lost by the donor. Consequently much of 
the donor's mass will be ejected from the binary system so 
that, in the described model at least, Case C mass trans- 
fer will be non-conservative. The mass lost from the sys- 
tem will still be relatively close by when the SN explodes 



© 2007 RAS, MNRAS 000, 000-000 



A deeper search for the progenitor of the Type Ic Supernova 2002ap 17 



since Case C mass transfer occurs late in the evolution of 
the donor star. Exactly what distance this material has dis- 
persed to depends on the velocity v with which it is ejected 
from the binary and the elapsed time t between the end of 
mass transfer and the explosion of the SN. This distance 
is given by D « 3.15 x 10 9 (u/100 kms^)(i/l yr) km (or D 
« 21(u/100kms~ 1 )(t/lyr) au). The velocity of the ejected 
material v will be similar to the orbital velocity of the ac- 
creting star, of the order of 10 - 100 kms -1 . In Section [5.21 
we used the radio observations of the SN to calculate the 
minimum radial distance out to which no significant density 
of CSM was encountered by the SN ejecta: D m m = 11,000 
to 25,000 au. Assuming v — 100 kms -1 requires that mass 
transfer ended at least t m i n = 500 to 1200 yrs prior to the 
Since this time period is rather short we cannot rule 
out all Case C mass transfer scenarios. Binaries with or- 
bital radii approaching the lower limit of what constitutes a 
Case C system might begin to interact early enough so as to 
satisfy the above requirement. The low extinction measured 
towards the SN (see Section |4| probably suggests that the 
10 to 15M© of material lost prior to explosion had dispersed 
to much larger distances than the lower limit we estimate 
from the radio and X-ray data. This would tend to rule out 
all Case C binaries. One can still, however, invoke a non- 
spherical geometry (e.g. the rings seen around SN 1987A) 
that would allow this material to be much closer without 
causing extinction of the SN provided the line-of-sight to 
the observer is not obscured. 

Case B mass transfer is where the donor star begins to 
transfer mass when it first becomes a giant at the end of core 
hydrogen burning, but before the ignition of central helium 
burning. Orbital periods of systems undergoing this type 
of interaction are intermediate to those of Cases A and C. 
As with Case C, the mass loss will occur on the relatively 
short thermal or dynamical timescales, making it difficult 
for a main sequence companion to accrete all this material. 
Case B mass transfer will happen early enough in the evolu- 
tion of the donor star to allow sufficient time for the ejected 
material to disperse prior to the SN explosion. This earlier 
stripping of its hydrogen envelope will also afford the donor 
star more time to shed its helium envelope through strong 
stellar winds and become the low mass WC star predicted 
to be the progenitor of this SN. Based on this argument we 
suggest that the mass transfer process in our assumed bi- 
nary system is more likely to have been Case B than Case 
C. We conclude that in all possible binary systems for the 
progenitor of SN 2002ap where the companion star was ini- 
tially less massive than the progenitor, the companion must 
have been initially ^ 20M© and non-conservative Case B 
mass transfer is most likely to have occurred. 

To test this idea we created binary models using the 
STARS stellar evolution code (Eldridge, Izzard & Tout in 
prep). We note that the modelling of binary stars is ex- 
tremely uncertain (e.g. Hurley et al. 2002b) . In all the mod- 
els the initial orbital periods were such that Case B mass 



° The model assumed here is rather simple. In reality the strong 
stellar wind of the newly formed W-R star would collide with the 
previously ejected material causing it to accelerate (Eldridge et 
al. 2006 and references therein). This would have the effect of 
reducing f m i n above. 



transfer would occur. One of the models, incorporating a 
2OM0 primary and a 14Mq secondary, is shown as an ex- 
ample in Figure [TJJ compared to the relevant luminosity lim- 
its. The evolution of both models continues up to the end 
of core carbon burning in the primary. At this point, the 
time until the primary goes SN is so short that both stars 
appear essentially as they will when it explodes. Both stars 
must fall below the luminosity limits at the time of the SN 
explosion in order to satisfy observations. Being of lower 
mass, the secondary (the companion) is still on the main se- 
quence as the primary expands and begins to transfer mass. 
However in this model, and in fact in all of the Case B bi- 
nary models created, it was found that mass loss through 
Roche lobe overflow could not remove mass quickly enough 
to halt the expansion of the primary star, with the result 
that it expanded to engulf its companion. Such a scenario is 
referred to as common envelope evolution (CEE). In this ex- 
ample CEE begins with the dense core of the primary and 
the main sequence secondary orbiting about the centre of 
mass of the system within the inflated envelope of the pri- 
mary. Orbital energy is transferred to the common envelope 
through an unknown physical process, leading to a reduction 
in the orbital period and radius, and to the ejection of the 
common envelope. The timescale for CEE is short and there- 
fore very little mass is accreted by the secondary during this 
time. Provided the binary orbit does not shrink so much as 
to cause a mergeiQ, having been stripped of its H-rich enve- 
lope the primary star continues to lose mass through strong 
stellar winds, eventually shedding most of its He-rich layer 
also. The model shown in Figure [TT] produces a 5.4M© C+O 
W-R from the initially 20M Q star. 

In the case of an assumed neutron star or black hole 
companion we derive broadly the same conclusions about 
mass transfer processes as we do for a main sequence sec- 
ondary. There are however some notable differences. Case 
A mass transfer would again most likely lead to the merger 
of the binary components. However the merger of a main 
sequence star with a neutron star or black hole would ei- 
ther completely disrupt the progenitor, or possibly result in 
the formation of a Thorne-Zytkow object (TZO) (Thorne 
& Zytkow 1977). Fryer & Woosley (1998) suggest that the 
merger of a black hole or neutron star with the He core of 
a red giant could produce a long complex GRB with no as- 
sociated supernova. This may occur in Case B systems with 
initially small orbital separations. In general both Case B 
and Case C mass transfer processes are possible, either with 
or without a subsequent CEE phase, although we again sug- 
gest that Case B is more likely since it results in the earlier 
removal of the progenitor's hydrogen envelope and a much 
longer period of W-R evolution prior to explosion. We are 
of course unable to place mass limits on a neutron star or 
black hole companion using our luminosity limits. 

Note that in all of the above we have employed binary 
interaction to remove only the hyrogen envelope of the pro- 
genitor star. The subsequent removal of the He-rich layer 
has occurred via strong radiatively driven winds. Nomoto et 
al. (2001) point out that under the right conditions a second 



7 Nomoto et al. (2001) propose that some binary mergers of this 
sort may produce rapidly rotating C+O stars, which could ex- 
plode as Type Ic hypernovae. 
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Figure 11. Binary model of a 2OA/0 primary and a IAMq secondary created using the Cambridge STARS stellar evolution code and the 
binary-evolution algorithm of Hurley et al. (2002b). The luminosity limits plotted are derived from the pre-explosion observations in the 
W-R domain, and from post-explosion observations for O to M-type supergiants. The initial separation is such that the primary fills its 
Roche lobe as it expands to become a giant after core hydrogen burning, and so Case B mass transfer ensues. However, mass loss through 
Roche lobe overflow does not remove mass quickly enough to halt the expansion of the primary star, with the result that the primary 
star expands to engulf its companion. The system has entered a common envelope evolution (CEE) phase. Orbital energy is transferred 
to the common envelope through an unknown physical process, leading to a reduction in the binary period and separation, and to the 
ejection of the common envelope. The evolution of both models is halted at the end of core carbon burning in the primary, essentially just 
prior to the SN explosion. The numbers displayed beside the primary and secondary tracks show initial and final masses. The primary 
ends its life as a 5AMq WC star. The secondary accretes very little of the ^15Mq of material lost by the primary, increasing in mass 
by just 0.06M Q . 



mass transfer event may occur, which removes the He layer. 
This however requires a large degree of fine-tuning of the 
binary system. During the first CEE phase the binary orbit 
must be reduced enough to allow interaction between the 
newly formed He-star and its companion, but not so much 
as to cause the two to merge. This second mass transfer is 
more likely to occur for lower-mass He-stars since they at- 
tain larger stellar radii (Habets 1986). Nomoto et al. (1994) 
exclusively use such double mass transfer to model the evo- 
lution of the progenitor of the Type Ic SN 19941; a 2.1M Q 
C+O star formed from a 4Mq He-star. The significantly 
larger estimates for the final progenitor mass of SN 2002ap 
and the He-core mass from which it formed (see Section l5,2p 
makes the double mass transfer scheme less likely in this 
case. 



6 CONCLUSIONS 

We have used a unique combination of deep, high qual- 
ity CFHT pre-explosion images of the site of SN 2002ap 



and follow up HST and WHT images of the SN itself to 
place the most restrictive limits to date on the luminosity 
and mass of the progenitor of a Type Ic SN. Theory pre- 
dicts that the progenitors of such SNe are W-R stars which 
have lost all of their H and most of their He-rich envelopes. 
The archival observations rule out as viable progenitors all 
evolved hydrogen-rich stars with initial masses greater than 
7-8Mq, the lower mass limit for stars to undergo core col- 
lapse. This is entirely consistent with the observed absence of 
hydrogen in the spectra of this SN and with the subsequent 
prediction that its progenitor was a W-R star. The magni- 
tude limits do not allow us to distinguish between WN, WC 
and WO stars, as examples of all these W-R types have been 
observed with magnitudes that fall below our detection lim- 
its. Since the SN was deficient of helium it is likely to have 
been a WC or WO star. Using the W-R Mass-Luminosity 
relationship we calculated the upper limit for the final mass 
of a W-R progenitor to be < I2I3MQ. Comparing the W- 
R luminosity limit with Geneva and STARS stellar models 
for single stars at metallicity Z=0.008 we found no viable 
single star progenitors when standard mass loss rates were 
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considered, and only initially very massive (> 30 — 40Mq) 
candidate progenitors when mass loss rates were doubled. 
The lack of significant extinction towards the SN, along with 
weak detections of the SN at X-ray and radio wavelengths, 
suggests a relatively low mass of circumstellar material sur- 
rounding the progenitor star. This would tend to rule out 
very high mass progenitors. We conclude that any single star 
progenitor must have experienced at least twice the stan- 
dard mass- loss rates, been initially > 30Mq and exploded 
as a W-R star of mass 1O-12M0. 

The most likely alternative to the single star models is 
a binary evolution scenario, in which the progenitor star is 
stripped of most of its mass by a companion, becoming a low 
mass, low luminosity W-R star prior to exploding as a super- 
nova. This idea is consistent with Mazzali et al. (2002) who, 
based on modelling of the supernova lightcurve and spectra, 
suggest a binary scenario for the progenitor of SN 2002ap in 
which a star of initial mass 2O-25A/0 becomes a 5Mq C+O 
star prior to explosion. (Note that we revise this initial mass 
to the lower value of 15-20Mq in Section \5. 21 1 The luminos- 
ity limits from both the pre- and post-explosion observa- 
tions, along with other observational constraints, allowed us 
to place limits on the properties of any possible binary com- 
panion, and to infer which mass transfer processes could or 
could not have occurred. We conclude that, if the system 
was indeed an interacting binary, the companion star was 
either a main sequence star of ^ 2OA/0, or a black hole or 
neutron star formed when an initially more massive compan- 
ion exploded as a supernova. In either case, mass transfer 
most likely began as non-conservative Case B, which led im- 
mediately to common envelope evolution (CEE). 
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